Symmetric and asymmetric states in a mesoscopic superconducting wire in the 

voltage driven regime 
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The response of a mesoscopic homogeneous superconducting wire, connected with bulk normal 
metal reservoirs, is investigated theoretically as function of the applied voltage. The finite relaxation 
length of the nonequilibrium quasiparticle distribution function Le is included where we assumed 
that our wire is in the dirty limit. We found that both symmetric and asymmetric states can 
exist which are characterized by a stationary symmetric and asymmetric distribution of the order 
parameter with respect to the center of the wire. Current voltage characteristics of the wire with 
length L > Le being in the symmetric state show a pronounced S-behavior. The asymmetric state 
may exist only for voltages larger than some critical value and coexist with the symmetric state 
in a finite voltage interval. For wires with L ~ Le the asymmetric state survives up to higher 
values of the voltage than the symmetric one and may exist both in the voltage and the current 
driven regimes. We propose an experiment to observe reversible switching between those stationary 
symmetric and asymmetric states. 

PACS numbers: 74.25.Op, 74.20.De, 73.23.-b 



I. INTRODUCTION 



Modifications of the quasiparticle distribution func- 
tion f(E) in a superconductor will drastically influence 
its properties. Injection(extraction) of the quasiparticles 
via a tunnel junction suppresses(enhances) the supercon- 
ducting gap [l|, 0] and creates a charge imbalance in the 
sample [3j; external electro- magnetic radiation excites 
quasiparticles to higher energy levels and can consider- 
ably increase the superconducting gap and the critical 
current 0, Q; moving vortices will change f(E) resulting 
in the dependence of the vortex viscosity on the vortex 
velocity @; fast oscillations of the superconducting gap 
in superconducting bridges or in phase slip centers leads 
to a dynamical nonequilibrium distribution f(E) and re- 
sults in strongly nonlinear current-voltage characteristics 

0- 

Recently, the effect of nonequilibrium f(E) induced by 
an applied voltage was studied in a dirty superconducting 
wire (whose mean path length I is smaller than the coher- 
ence length at zero temperature £o) that was connected 
to large normal metal reservoirs with no tunnel barriers 
Q. Theoretically it was found that there is a critical 
voltage V c above which the superconductor exhibits a 
first order transition to the normal state. This effect is 
connected with a strong modification of the quasiparti- 
cle distribution function f(E) and especially its odd part 
f L (E) = f(-E) - f{E), by the applied voltage. The 
applied voltage leads to the excitation of quasiparticles 
to energy levels near the energy gap resulting in a strong 
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modification (i.e. diminishing) of the order parameter. 

In the present work we study the same system as 
Ref. Q but include the finite relaxation time te of the 
nonequilibrium quasiparticle distribution function f(E) 
due to inelastic electron-phonon interaction which was 
neglected in Ref. Q. We show, that the finite length 
Le = V t eD (D is the diffusion constant) weakens the 
effect of the applied voltage on the superconducting prop- 
erties and changes qualitatively the shape of the IV char- 
acteristic (it becomes S-shaped) already for wires with 
length L of several Le- Furthermore we find that in our 
homogeneous wire there are stationary states with both 
symmetric and asymmetric distributions of the order pa- 
rameter (with respect to the center of the wire). The 
asymmetric state exists only at voltage larger than some 
critical value and can coexist with the stationary sym- 
metric states in a finite voltage interval. For wires with 
length L ^ Le it survives up to larger applied voltages 
than the symmetric one and at specific conditions even 
for larger currents. 

The paper is organized as follows. In section II we 
present our theoretical model. In sections III and IV 
we discuss the existence of the symmetric and asymmet- 
ric states, respectively. Finally, in section V we present 
our conclusions and propose an experimental setup to 
observe the reversible transitions between the predicted 
symmetric and asymmetric states. 



II. MODEL 

We used a quasi-classical approach to calculate the 
nonequilibrium properties of a superconductor in the 
dirty limit and restrict ourselves to temperatures near 
T c . This allows us to use the Usadel equation 0] for 
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the normal a(E) — cosO and anomalous (3(E) = sinG 
quasi-classical Green functions in relatively simple form 

^-f + ((2i J B--^-)-(V^) 2 cose)sine + 2|A|cose = 0. 

&X E 

(i) 

In the same limit the diffusive type equation for the 
space dependence of the transverse (even in energy) 
f T (E) = 1 - f(E) - f(—E) and longitudinal (odd in 
energy) fh{E) = f(—E) — f(E) parts of the quasiparti- 
cle distribution function 2f(E) = 1 — fh{E) — fr(E) are 
given by 

V((^ 2 - R 2 2 )Vf L ) + 2N 2 R 2 VipVf T 

-w(/i-/E)=0, (2a) 

V((JV? + A/f )V/ T ) + 2N 2 R 2 V V Vf L - 
N elf® 

ji(f T -^)-2N 2 \A\f T = 0. (2b) 

E 

Here <p is the phase of the order parameter A = |A|e J<?i , (j> 
is an electrostatic potential, Ni(E) + iRi(E) — cosQ(E), 
N 2 (E) + iR 2 (E) = sine{E) and /£(£) = tanh(S/2T) 
is the odd part of the equilibrium Fermi-Dirac distribu- 
tion function of the quasiparticles. The dimensionless 
length Le = V Dte /£o = ^teAq/H defines the range 
over which the nonequilibrium distribution of the quasi- 
particles may exist in the sample. 

Within the same approach we have the rather simple 
self-consistent stationary equation for the order parame- 
ter A 

ai^ + (l-a 2 |A| 2 + * 1 +z* 2 )A = 0, (3) 

which is analogous to the Ginzburg-Landau equation [? 

but with the additional terms \&i = f °° R 2 (fL — 
f L )dE/\A\ and * 2 = / °° N 2 f T dE/\ A\. Because we are 
interested in a stationary solution we removed all time- 
dependent terms in Ecl (2-3) which are present in the 

original equations. [ToL llll [l2| 

In Eqs. (1-3) the order parameter A is scaled by the 
zero-temperature value of the order parameter Ao ~ 
1.76fcsT c (in weak coupling limit), distance is in units 
of the zero temperature coherence length £ = \/hDjAo 
and temperature in units of the critical temperature T c . 
Because of this choice of scaling the numerical coefficients 
in Eq. (3) are ai ~ 0.69 and a 2 ~ 0.33. The current is 



scaled in units of jo — Ao/(£oPne) and the electrostatic 
potential is in units of </> = A /e (p n is the normal state 
resistivity and e is the electric charge). 

The deviation of /l(E) from its equilibrium value /£ 
may considerably influence the value of the order param- 
eter through the term ^i. The nonzero $t{E) is mainly 
responsible for the appearance of the charge imbalance in 
the superconductor and the conversion of the normal cur- 
rent to the superconducting one and vice versus. The lat- 
ter one occurs due to the Andreev reflection process (term 
27V2|A|/t in Eq. (2b)) or/and due to inelastic electron- 
phonon interaction (term Ni(/t — 4>df° /dE)/L 2 E in Eq. 
(2b)). 

The system of Eqs. (1-3) was numerically solved using 
an iterative procedure. First Eq. (3) was solved by the 
Euler method (we add a term dA/dt in the right hand 
side of Eq. (3)) until |A| became time independent. Than 
Eqs. (1-2) were solved from which we obtained new po- 
tentials and ^2 which were then inserted in Eq. (3). 
This iterations process was continued until convergence 
was reached. In some voltage interval no convergence was 
reached and we identify it as the absence of stationary 
solutions to Eqs. (1-3). Thus our numerical procedure 
automatically checks the stability of the found stationary 
symmetric and asymmetric states. 

The numerical solution of Eqs. (1-3) converges much 
faster than the equivalent but more complicated Eqs. 
(5,7,8,9) in Ref. [a]. This allowed us to consider on 
the one hand a rather large interval of parameters but 
on the other hand it restricts the validity of the ob- 
tained results to the temperature interval 0.9 < T/T c < 
1 which roughly corresponds to the validity condition 
A(T)/k B T c < 1 of Eqs. (1-3) 

Knowing the solution of Eqs. (1-3) the current in the 
system can be found using the following equation 

j = 2 ai |A| 2 V<^ 

+ / ((iV 2 + Nl)Wf T + 2N 2 R 2 f L V V ) dE. (4) 
Jo 

The first term in Eq. (4) may be identified as a super- 
conducting current and the second one as a normal one. 
Assuming that the effect of the free charges is negligi- 
ble in the superconductors the electrostatic potential is 
determined by the following expression 

<$>= \ NxfrdE. (5) 
Jo 

We used the following boundary conditions for the system 
of Eqs. (1-3) 
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Q(±L/2) = A(±L/2) = 0. (6a) 
M ±m - i (tab (*±^M) + ...nh (^™)) , (6b) 

M ± L/2) . I (tanh (*±^) - tab (^™)) , (6c) 



which models the situation where there is a direct elec- 
trical contact of our superconducting wire of length L 
with large normal metal reservoirs at an applied voltage 
V{±L/2). 



III. SYMMETRIC STATES 

In Fig. 1 we present our calculated current-voltage 
(IV) characteristics of the superconducting wire with 
length L = 100, temperature T = 0.98 and different Le- 
In this section we consider the case of symmetric applied 
voltage V(±L/2) =tV. 

When L E < £(T) < L the f L (E) decays on the scale 
of the coherence length and practically does not influ- 
ence the value of the order parameter. When the current 
density in the wire exceeds the critical value (marked by 
Vq in Fig. 1), the N-S borders (formed near the ends 
of the wire due to the boundary condition Eq. (6c)) be- 
come unstable and start to move to the center of the wire. 
The superconducting region decreases and it effectively 
increases the resistance of the whole sample leading to 
a decrease of the current up to its critical value. This 
process results in the "vertical line" -like region in the IV 
characteristics of the wire. The longer the sample, the 
wider will be this "vertical line" -like region. Finally, at 
the critical voltage V c the region where | A| ^ is reduced 
to about £(T) and a first order phase transition to the 
normal state occurs. 

In this limit we may use the so called 'local equilib- 
rium model' which leads to the extended time-dependent 
Ginzburg-Landau (TDGL) equations [HI, [Hj]. A simple 
analysis based on those equations [l3| shows an increas- 
ing critical current with a stable N-S border when in- 
creasing the parameter 7 = 2teA /H — 2L 2 E . This is the 
reason for the increase of the maximal possible current 
in Fig. 1(a) for small values of Le- A direct comparison 
with results calculated within the extended TDGL equa- 
tions (see Fig. 2) show qualitatively the same behavior. 

In the opposite limit Le ~ L 3> £(T) the nonzero 
/l (E) considerably affects the value of the order param- 
eter everywhere in the sample. In Ref. @ it was shown 
for T = and Le = 00, that when the applied voltage is 
of about Ao the superconducting order parameter van- 
ishes everywhere in the sample and it leads to a first order 
phase transition to the normal state. Actually, qualita- 
tively the same is true for arbitrary temperature: when 
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FIG. 1: Current-voltage characteristics of the superconduct- 
ing wire with length L — 100, temperature T = 0.98 and dif- 
ferent relaxation lengths Le- Figure (a) corresponds to the 
case Le < £(T) ~ 7 and figure (b) to the case Le 3> £(T). 
Solid line shows the IV characteristic in the normal state. 
Dashed lines in both (a) and (b) mark the voltage interval 
where there are no stationary solutions to Eqs. (1-3). 



the voltage reaches approximately A(T) the order pa- 
rameter starts to decrease rapidly. In Fig. 3 we show the 
dependence of ]A|(V) Q for a spatially homogeneous 
distribution of A and /l(E) as found from a numerical 
solution of the model equation for the order parameter 
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FIG. 2: Comparison of the results obtained through solution 
of Eqs. (1-3) with results of a numerical solution of the ex- 
tended TDGL equations. Up to Ie ~ fCO t ne extended 
TDGL equations gives qualitatively the same results. Dashed 
lines mark the voltage interval where no stationary solutions 
of Eqs. (1-3) and the TDGL equations can be found. 
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FIG. 3: Dependence of the order parameter on the parameter 
V as obtained from Eq. (7) for different temperatures. 



with 8 = N{0)v/A Q = 0.3 (7Y(0) is the density of states 
of electrons at Fermi surface and v is a coupling constant) 
and where instead of the equilibrium weight factor due to 
the quasiparticles tanh(i?/2T) [? ], we use the nonequi- 
librium one for the case of half voltage drop V. 

At low temperatures the order parameter practically 
does not vary with V up to some critical voltage and it 
results in the absence of that parts in the IV character- 
istics with negative differential resistance (see Ref. 
Contrary, at higher temperatures, |A| may vary strongly 
with V (see Fig. 3) and it leads to a decay of the full 
current in the sample because effectively the resistance of 
the wire grows with increasing the voltage. It results in 
the appearance of a pronounced part with negative differ- 
ential resistance dV/dl in the IV characteristic (see Fig. 
2 in Ref. @ and Fig. 1(b) for L% = 2 12 and L% = 2 10 ) 



for V > Vq. No stationary solutions to Eq. (1-3) were 
found for samples with length larger than approximately 
10£(T) at V > Vq and temperatures not far from T c (see 
illustration in Fig. 4 for T = 0.98). This is connected 
with the fact that the superconducting current reaches a 
value close to the depairing current density in long sam- 
ples but for shorter ones the negative second derivative 
d 2 \A\/d 2 x starts to play an essential role in Eq. (3) and 
stabilizes the stationary state with nonzero |A|. 
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FIG. 4: Current-voltage characteristics of the superconduct- 
ing wire with relaxation length Le = 64, temperature T = 
0.98 and different wire lengths L (dashed line marks the volt- 
age interval with no stationary solutions of Eqs. (1-3)). It 
demonstrates the absence of the stationary states at V ~ Vb 
with increasing length of the wire. 

For intermediate wire lengths L ^> Le ^ £(T) the 
situation is more complicated. At V — Vq due to the 
suppression of the |A| near the edges by fb{E) the order 
parameter decays and the current decreases with increas- 
ing voltage. In Fig. 5 we show the spatial distribution 
of the order parameter for different values of the voltage 
for L E = 64 and L E = 16 (L=100, T=0.98). Due to the 
decay of Jl on the scale Le and the term in Eq. (3), 
the order parameter also starts to vary on a distance of 
about Le- That's why the order parameter is finite in 
a wider region of the wire with L ~ Le than for wires 
with L ^ Le even at V ~ V c . The suppression of |A| 
near the ends of the short wire (with respect to Le) is 
mainly connected with the relaxation of /l(E) due to 
the coupling with the transverse mode in Eq. (2a) (see 
discussion below). When the region of the suppres- 
sion of |A| increases up to about Le the effect of the 
nonequlibrium /l(E) becomes less pronounced, |A| de- 
creases much slower with increasing voltage and the cur- 
rent starts to increase again. At some moment it reaches 
a value close to the depairing current density in the cen- 
ter of the wire and the state with a nonequlibrium and 
a stationary distribution of the order parameter becomes 
unstable [Tg|. Phase slip centers should appear in the 
wire. This is observed as the absence of a stationary 
solution to Eqs. (1-3) in some range of applied voltage 
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FIG. 5: Distribution of the order parameter in a supercon- 
ducting wire of length L = 100 being in the symmetric su- 
perconducting state at different values of the applied voltage 
and relaxation lengths Le- 
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FIG. 6: Current-voltage characteristics of the superconduct- 
ing wires of different lengths and electron-phonon relaxation 
time (solid curves correspond to the normal state of the wires 
with lengths L = 100 and L = 200). It shows, that for a fixed 
temperature (not too close to T c in order to satisfy the con- 
dition Le 3> £(T)) and relaxation length Le it is possible to 
find the length of the sample for which the IV characteristic 
will have a pronounced S-behavior. 



(marked dashed lines in Fig. 1(b)). With further increas- 
ing voltage the superconducting region decreases and we 
find again a stationary solution. At this voltage, the size 
of the superconducting region is too small to 'carry' a 
phase slip state and is close to the critical size 10£(T) 
found above for the stability of the stationary state in 
the wire with Le ~ L. The above mechanism leads to 
a S-shape for the IV characteristics for the intermediate 
case L ^> Le £,{T). This feature is more visible if we 
consider low temperatures and study samples of different 
lengths (see Fig. 6). 

The relaxation of fL(E) and fr(E) to equilib- 



rium occurs not only via the inelastic electron- 
phonon interaction but also due to the coupling terms 
2LW2-R2V</?V/t,l in Eqs. (2a, b). They are nonzero in 
the region where the order parameter changes near the 
N-S border and where V/t,l 7^ 0. This mechanism ef- 
fectively relaxes the nonequilibrium Jl,t{E) at energies 
close to |A| (where the product N2R2 7^ 0). In Fig. 7 we 
present the spatial dependence of Jl(E,x) — f L (E) and 
fr(E, x) for different values of the energy for a wire with 
length L=100, temperature T=0.98, relaxation length 
Le = 64, applied voltage V = 0.08 and with zero and 
nonzero coupling terms 2£W 2 i?2V<y9V/T,L in Eq. (2). 
We notice that for energies close to |A| the coupling of 
the longitudinal mode with the transverse one leads even 
to a change in sign of fL(E) in the central part of the 
sample and hence to a smaller value of the term \E'i in Eq. 
(3) . The change of the relaxation length of the transverse 
mode fx is not so pronounced, but is still noticeable and 
justifies the faster decay due to coupling with the longi- 
tudinal mode. 
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at different energies for zero and nonzero coupling terms 
2N2R2 V(^V/l,t in Eqs. (2a, b). Parameters of the wire: 
L=100, L| = 2 9 , T = 0.98, V = 0.08, A{x = 0, V = 0.08) ~ 
0.216. 



6 



That effect provides the dependence of the shape of 
the current-voltage characteristic of the wire on temper- 
ature (see Fig. 8) for fixed length of the sample and 
strength of the electron-phonon interaction even when 
L ^> Le > £,(T). Above discussed additional mech- 
anism of relaxation effectively shortens the relaxation 
length of /l.t(E) for V > Vq and decreases the effect 
of nonequilibrium when Le approaches £(T). It moves 
the wire from the limit L ^> Le S> £,(T) to the limit 
L^> L E ~ at some temperature close to T c . 
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FIG. 8: Current-voltage characteristics of the superconduct- 
ing wire of fixed length and strength of the electron-phonon 
interaction for different temperatures. It shows the change of 
the shape of the IV characteristic when Le approaches f(T). 



IV. ASYMMETRIC STATES 

To understand the origin of the asymmetric states in 
a superconducting wire let us consider as a example the 
mesoscopic normal metal wire half of which has a resis- 
tance larger than the other half. Applying a voltage to 
such a sample the voltage drop in the high resistance 
part will be larger than in the low resistance part (for 
definition we suppose that Ri e ft > Rright) in order to 
satisfy the continuity of the current along the wire, i.e. 
I(x) = const. It means that if we put V=0 in the cen- 
ter of our inhomogeneous wire we have \V(— L/2)\ > 
\V(L/2)\. Hence the deviation of the quasiparticle dis- 
tribution functions from equilibrium in the left and right 
normal reservoirs (f(E,x±L/2) — l/(exp(.E/T) + l)) will 
depend not only on the value of the applied voltage but 
also on the properties of the mesoscopic wire. 

To find corresponding superconducting state in our ho- 
mogeneous system we use the asymmetric boundary con- 
ditions V(±L/2) = TV + dV. The addition dV is deter- 
mined from the condition of constant current along the 
wire. It increases from zero (at V ~ Vo) to a finite value 
(at V£ - critical voltage of the asymmetric state). No 
asymmetric states were found for V < Vq. In Fig. 9 we 
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(empty symbols) and asymmetric (filled symbols) states. 
Solid curve shows the normal state. Dashed curves mark re- 
gions of absence of the stationary symmetric states. The sta- 
tionary asymmetric state does not exist for V < Vo - voltage 
where negative-differential resistance appears for symmetric 
state. 



show typical IV characteristics and in Fig. 10 we present 
the corresponding distribution of the order parameter for 
symmetric and asymmetric states at the same values of 
the voltage. It is obvious from the above discussion that 
the order parameter is more suppressed on the side of 
the wire with maximal absolute value of the voltage. For 
example if we change the sign of dV the order parameter 
distribution in Fig. 10 will be symmetrically reflected 
with respect to x=0. 

It is interesting to note that such an asymmetric 
state may exist up to larger values of the voltage when 
L > Le S> £, and even for larger values of the current, 
than the symmetric state (latter property is realized for 
lengths L ~ Le for which the IV characteristic changes 
slope from a negative to a positive one for the symmetric 
state when V > Vq). Another interesting property is the 
existence of the stationary asymmetric state for voltages 
when the stationary symmetric state is absent (see Fig. 
9). 

The range of the existence of the stationary asymmet- 
ric state is rather small in comparison with the symmetric 
state when the length of the wire is much larger than Le- 
It exists only at voltages near Vo (see Figs. 9 and 11 for 
small ratios Le/L). 

We did not find stationary asymmetric states for wires 
with Le ~ £(T). For small values Le <C the asym- 
metric state is nearly the same as the symmetric one 
but shifted with respect to the point x — and the IV 
characteristics are found to be nearly the same for both 
symmetric and asymmetric states. This is explained by 
the negligible effect of /l(-E) on the value of the order 
parameter in this limit. 
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FIG. 10: Order parameter distribution for asymmetric and 
symmetric states. 
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FIG. 11: Current- voltage characteristics for symmetric and 
asymmetric states for wires with different Le- 



V. CONCLUSION 

The response of a dirty superconducting wire that is 
connected to normal metal banks, to an applied volt- 
age, strongly depends on the ratio between the length 
of the sample L, the coherence length £ and relaxation 
length Le of the nonequilibrium quasiparticle distribu- 
tion function f(E). The latter one, was determined in 
our calculations by the inelastic electron-phonon inter- 
action Le = V Dte with an energy independent char- 
acteristic time te- Applied voltage suppresses the su- 
perconducting order parameter both via the creation of 
the superconducting current and the modification of the 
quasiparticle distribution function. First mechanism is 
prevalent when Le < £(T) and the second one is domi- 
nant in the opposite limit Le ^> £,(T) for samples with 
L ~ Le- For wires with L ^> Le S> the transition 
to the normal state is mainly induced by the supercon- 



ducting current. 

We found that for the following two cases L ^> £(T) > 
Le and L > Le S> there are three stationary so- 

lutions to Eqs. (1-3) at given value of the applied volt- 
age difference 2V > 2Vq - one symmetric state and two 
asymmetric states (which are symmetrically reflected to 
the center of the wire) which are characterized by sym- 
metric and asymmetric distribution of the order param- 
eter, respectively. The degeneracy is most pronounced 
in the limit L > Le S> £(T) because the order param- 
eter may vary on the distance of order Le S> £ due to 
the long relaxation length of the odd (in energy) part 
of the nonequilibrium quasiparticle distribution function 
and it provides the basis for the appearance of new ef- 
fects. For example it leads to an S-behavior of the IV 
characteristic of the wire being in the symmetric state. 
In the same limit the stationary asymmetric state may 
exist even when the symmetric one exist only as a time- 
dependent one (phase slip state) (see Fig. 9). Moreover, 
for a specific temperature and length Le S> £(T) there 
exist length of the sample for which a stationary asym- 
metric state does exist up to a larger value of the current 
than the stationary symmetric one (see Fig. 9 and com- 
pare for L% = 2 s ). 

We should mention here, that in the extended TDGL 
equation there is no direct effect of the voltage on the 
value of the order parameter. This is because of the 
"local equilibrium approximation" (valid in the limit 
L E < £(T)) [HG3 where it was assumed that f L - f° 
is only proportional to the time variation of the value 
of the order parameter d\A\/dt and hence is zero in 
the stationary case. But it is interesting to note that 
the asymmetric state could also be obtained by this 
equation. For example the asymmetric state exists for 
7 = 32, L = 100, T = 0.98 (see inset in Fig. la) at 
0.217 < V < 0.242. But we did not find in this model 
the wire parameters for which asymmetric and symmet- 
ric states could exist at the same value of the voltage. 

Experimentally the asymmetric state may be realized 
by adding additional contacts to the superconducting 
wire with control current close to the ends of the wire 
(see Fig. 10). It is better to contact the wire with the 
reservoirs made from the same material and to apply a 
strong enough magnetic field to suppress the supercon- 
ductivity in the reservoirs. This procedure will provide 
good contact with normal reservoirs. Applying strong 
enough current to only one of the current contacts we lo- 
cally destroy superconductivity and force an asymmetric 
distribution of the order parameter. After switching off 
the control current the asymmetric distribution should 
be stable (at proper choice of the working point at the 
IV characteristic). To come back to the symmetric state 
we should apply a control current at the both current 
contacts and it will force the recovering the symmetric 
state in the wire. We will need a series of tunnel junc- 
tions along the wire to measure locally the strength of 
the order parameter. Other alternative is to use super- 
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Tunnel junctions 
or superconducting probes 



FIG. 12: Proposed experimental setup to measure transi- 
tions between symmetric and asymmetric states and current- 
voltage characteristic of both states. The curves with empty 
and filled symbols show the order parameter distribution for 
a wire with parameters from Fig. 10 and V=0.389. 



relatively large coherence lengths (£a;(Q) ~ 0.15/irn, 
£zn(0) ~ 0.25/im ) and Le(AI) ~ 9/im, Le{Zti) ~ 35/im 
(L E (Al) ~ 60, L E (Zn) ~ 140). The expected tempera- 
ture range for the validity of our results is 0.9 < T/T c < 
1. The most strict restriction to observe these effects is 
that there should be only a small heating of the sample. 
Heating will initiate the transition of the sample to the 
normal state and hide the main effects. 

For other low-temperature superconductors (Nb, Pb, 
In, Sn) we have Le < 10 (we used data for te at T ~ T c 
from Ref. [l6j]). It means that at T > 0.9 the coher- 
ence length is comparable with L e ■ For these conditions 
we did not find asymmetric states. Moreover the effect 
of nonequilibrium in the odd part of the quasiparticle 
distribution function that is connected with the applied 
voltage is rather small for these parameters and the ex- 
tended TDGL equations gives result which are qualita- 
tively similar to the one obtained here. 



conducting (SC) probes to distinguish the part of the 
superconductor with attached SC-probe in the normal 
or superconducting state. When in the normal state the 
superconducting probe measures the electrostatic poten- 
tial. 

Good candidates to observe these symmetric and 
asymmetric states is dirty aluminium and zinc with their 
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